#### Figure 5: Bivariate Support for Loyalty ####

# Libraries
# library(here)
# library(rio)
# library(tidyverse)
# library(ggthemes)
# library(lubridate)
# library(estimatr)

# data_pnas = import(here("Data","data_pnas.rds"))
# source(here("Code","Functions","funcs_theme.R"))
# source(here("Code","Functions","funcs_bvplot.R"))

fig_5 = plot_bv_maga(data_pnas, norm_loyaltyc)

print(fig_5)

ggsave(here("Plots","Main","figure_5.pdf"),
       fig_5,
       units = "in",
       width = 6, height = 6,
       dpi = 600)